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1.  Introduction 


The  point  spread  function  (PSF)  is  the  image  response  of  an  imaging  system  in  sensing  a  scene 
containing  a  single  point  scatterer.  As  an  extended  scatterer  can  be  approximated  as  a  collection 
of  point  scatterers  in  many  cases,  the  PSF  is  useful  for  characterizing  the  image  response  in 
generalized  sensing  scenarios.  The  PSF  then  provides  a  mathematical  framework  for  quantifying 
the  effects  of  the  sensing  geometry  and  associated  parameters — such  as  aperture  configuration 
(single-  vs.  multi-aperture),  array  aperture  size,  spectral  bandwidth,  spatial  and  spectral  sampling 
interval,  etc. — on  the  imaging  and  detection  performance.  Specifically,  the  PSF  offers  guidance 
on  1)  the  design  and  optimization  the  imaging  system,  2)  the  prediction  of  the  detectability  limit 
of  the  system,  and  3)  the  interpretation  of  the  scene  image  response. 

In  this  work,  the  PSF  is  studied  within  the  context  of  ultra-wideband  (UWB),  electromagnetic- 
oriented  standoff  target  localization.  The  twofold  objective  here  is  to  first  derive  efficient, 
simplified  closed-form  approximations  for  the  single-  and  multi-aperture  forward-looking  PSF 
for  on-surface  and  subsurface  scatterers  and  then,  using  those  approximations,  put  forth 
observations  useful  for  understanding  the  imaging  results.  Although  the  focus  of  the  study  is  on  a 
sensing  setup  that  emulates  the  operation  of  an  existing  forward-looking  radar  system,  it  is 
expected  the  derivations  featured  can  be  generalized  for  the  treatment  of  a  variety  of  standoff 
sensing  configurations. 


2.  Analytical  Formulations 


From  the  time-reversal-based  imaging  functional,  the  PSF  is  first  formulated  for  an  on-surface 
scatterer  and  then  the  derivation  is  extended  for  a  subsurface  scatterer. 


2.1  On-Surface  Scatterer 


The  sensing  geometry  treated  here  is  akin  to  the  one  studied  in  Reference  1 .  Consider  an  array  of 
vertically  polarized  point  radiators  and  a  near-ground  or  on-surface  point  scatterer  at 
r  =  (r  ,  9 s  ,(/>s),  approximating  the  Green’s  function  for  the  interactions  of  the  co-polarized  field 

component  as2 
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and  assuming  observation  with  an  array  of  isotropic  point  receivers,  the  PSF  for  a  single-aperture 
configuration  (positioned  above  the  ground  at  height  zT  =  zR)  can  be  constructed  from  the  time- 

3 

reversal-based  imaging  functional  as 
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Using  far-field  expressions  for  the  distance  vectors,  the  F (•)  terms  can  be  moved  outside  the 

summation  operation,  and  orienting  the  transmitting  and  receiving  arrays  parallel  to  the  v-axis 
( xT  =xR  =  0  ),  the  function  in  Eq.  2  becomes 
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where  8k  ,  SyR ,  and  SyT  are  the  A-space  sampling  interval,  the  receiver  spacing,  and  the 
transmitter  spacing,  respectively;  kc  is  the  center  wavenumber;  the  wavenumber  bandwidth  is 
then  A k  =  Nk8k ,  and  the  aperture  sizes  are  LR  =  NR8yR  and  Lr  =  NT8yT .  Applying  a 
stationary-phase-like  approach  to  evaluate  the  summation  in  A-space  and  noting  the  relation 
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the  amplitude  of  Eq.  3  approximately  simplifies  to 
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where  it  is  assumed  that  there  is  a  transmitter  at  each  end  of  the  aperture  LT .  The  first 
sin(-)  /  sin(-)  term  in  Eq.  5  dictates  the  behavior  of  the  PSF  in  range:  on  the  xv-plane 
( 6  =  0s  =  k  /  2 ),  it  can  be  shown  that  the  first  null  occurs  at  cl  2 A f  — which  is  the  resolution  in 
range,  8r ,  independent  of  r  and  (f>s .  The  periodicity  of  this  tenn  detennines  the  frequency 
sampling  interval  needed:  for  maximum  range  r  max ,  the  condition  Sf  <  c  /  2rs  max  must  be  met. 

Similarly,  the  cross-range  properties  of  the  PSF  are  controlled  by  the  last  two  terms  in  Eq.  5.  For 
the  current  two-transmitter  system  (and  LT  =  LR  ),  based  on  the  location  of  the  first  null  of  these 

two  terms,  it  can  be  readily  seen  that  the  angular  resolution  depends  on  the  angular  position  of 
the  scatterer  as 
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At  =  0 ,  =  Xc  /  2Zr  ;  the  angular  resolution  degrades  as  the  scatterer  moves  away  from  the 

broadside  direction  ( (f>s  =  0  )  of  the  aperture.  In  general,  the  spatial  sampling  should  be  less  than 
Xmin  /  2 ;  however,  for  an  UWB  system,  it  has  been  shown  that  high  cross-range  resolution  can  be 
achieved  with  coarser  sampling  than  the  stated  criterion.4 

For  transmitters  and  receivers  with  symmetric  patterns  about  the  array  axis  (or,  in  this  case,  the 
v-axis),  each  scatterer  would  give  rise  to  two  identical  image  responses — one  at  the  true  location 
of  the  scatterer  (at  <f>s )  and  an  ambiguous  component  at  the  location  that  is  mirrored  about  the 

array  axis  (at  7t-(j)s).  The  ambiguous  component  can  be  suppressed  through  the  use  of  the  multi¬ 
aperture  sensing  geometry — described  below — or  by  first  employing  transmitters  and  receivers 
with  forward  directivity  and  then  explicitly  incorporating  the  antenna  pattern  functions  into  the 
imaging  reconstruction.  Both  of  these  methods  would  reduce  the  mirrored  image  components 
originated  from  targets  positioned  in  the  backside  of  the  array,  which  can  mask  target  responses 
in  the  forward-looking  image. 


Following  a  line  of  derivations  similar  to  that  summarized  above,  the  approximate  PSF  for  the 
multi-aperture  configuration  with  interweaved  left  and  right  excitation  scheme  can  be  found  as 
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in  which  Lx  is  the  forward-traveling  distance  and  Sx  the  separation  between  successive  left  or 
right  apertures  ( L  «  Nx8x  ).  For  the  set  of  parameters  commonly  employed  in  forward-looking 
studies,  in  determining  the  range  and  cross-range  imaging  behaviors,  the  effect  of  the  terms 
(^(sin^cos^-sin^cos^M  and  r>xkc(sin<9cos0-sin6(cos^)/2  can  be  considered  to  be 
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negligible.  As  such,  the  frequency  sampling  criterion  and  range  resolution  are  the  same  as  those 
for  the  single-aperture  case.  The  cross-range  resolution  is  dependent  on  the  angular  position  of 
the  scatterer  as  before;  specifically,  again  on  thexy-plane,  S(f>MA  =  min  (<74,,  (5#  j ,  with 
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For  Lx  >  Lt  (e.g.,  Lx  =10  m  and  Lr  =2  m),  higher  angular  resolution  occurs  over  the  region 
away  from  the  broadside  direction  of  the  aperture:  at  #  =  0 ,  S(f>MA  =  Xc  /  2 LT  ,  and  the  highest 
resolution  is  reached  at  #  =  n  /  2 ,  with  8(f) MA  =  Xc  /  2Lx . 


The  resolution  in  9  of  both  the  single-  and  multi-aperture  cases  can  be  deduced  based  on  the  first 
null  of  the  first  sin(-)  /  sin(-)  term  in  Eqs.  5  and  7.  Upon  settingr  =  r. ,  it  is  found  that 

89 SA  =  80 ma  =  AT  /  2zr,  (9) 


where  AA  =  c  /  A/ .  The  above  result  stipulates  zr  is  relatively  large  (e.g.,  zT  =  2  m);  otherwise, 

the  other  tenns  within  the  PSF  would  also  need  to  be  taken  into  account.  Eq.  9  gives  the 
resolution  in  9  with  respect  to  the  coordinate  system  origin;  a  more  appropriate  definition  for  the 
resolution  in  elevation  can  be  derived  with  respect  to  the  aperture  center.  Designating  9'  as  the 
elevation  angle  with  respect  to  the  aperture  center,  for  a  multi-aperture  configuration,  with  the 
condition  that  0's  >  n  /  2  but  located  in  the  vicinity  of  the  xv-plane,  it  can  be  shown  that 
SO'MA  =  m\n(89'  ,89'  )  ,  in  which 
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The  above  supposes  the  terms  of  the  PSF  in  9'  corresponding  to  the  Lr  and  Lx  apertures  have 

definable  nulls;  if  only  one  of  these  nulls  is  definable,  then  that  null  is  used  to  detennine  the 
resolution.  Thus,  for  instance,  at  #  =  0 ,  89'MA  =  89[  ,  and  at  #  =  n  /  2 ,  80'MA  =  80]  .  These 

results  can  be  reinterpreted  for  the  single-aperture  case.  Specifically,  it  is  seen  that  at  #  =  0 ,  the 
single-aperture  case  does  not  have  resolution  in  O' ,  and  at  0  <  #  <  n  /  2  ,  S0'SA  =  89[t  — if  a  null 
is  definable. 


A  note  of  caution  needs  to  be  added  with  regard  to  the  suppression  of  the  imaging  sidelobes  in 
range.  In  the  setup  of  the  current  problem,  surface  clutter  close  to  the  transmitter  generates 
sidelobes  in  range  that  are  especially  detrimental  to  imaging  perfonnance.  This  is  due  to  the  fact 
that  scatterers  local  to  the  transmitter  not  only  produce  the  strongest  responses  but  also  create 
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sidelobes  that  spread  out  in  all  azimuth  directions.  Because  of  the  fast-decaying  nature  of  the 
backscattering  response  from  a  down-range  scatterer  (that  is,  asymptotically,  its  normalized 
backscattering  amplitude  is  of  the  form  cc  1  /  r2 ),  judicious  selection  of  both  5  f  and  a  frequency 
domain  window  is  necessary  to  ensure  the  sidelobes  of  the  near-range  surface  clutter  are  weaker 
than  the  response  from  the  down-range  scatterer.  It  is  seen  that  this  is  particularly  critical  for 
sensing  with  horizontal  polarization  since — as  evident  in  the  imaging  results1 — scatterers  have  a 
weaker  response  at  horizontal  polarization  than  at  vertical  polarization.  According  to  Eqs.  5  and 
7,  the  sidelobes  of  the  near-range  surface  clutter  should  fall  off  as  oc  1  /  r  ;  however,  their 
dependence  on  distance  in  the  imaging  algorithm  as  used  in  this  study  is  actually  oc  r  — as  a 
result  of  the  application  of  additional  amplitude  nonnalization  terms  needed  to  translate  the 
imaging  intensities  into  radar  cross  section  (RCS)  values.  In  view  of  the  above  observations,  a 
Blackman  window  with  a  =  0.08 5  is  chosen  and  imposed  on  the  frequency  domain  data.  A 
desirable  compromise  among  main  lobe  width,  sidelobe  amplitude,  and  sidelobe  decay  rate  is 
achieved  with  this  window.  Consistent  and  acceptable  results  are  obtained  even  though  the 
sidelobes  of  the  near-range  clutter — after  windowing — effectively  fall  off  only  as  cc  1  /  r  — a 
decay  not  as  fast  as  the  aforementioned  cc  1  /  r2  behavior  of  the  response  from  a  down-range 
scatterer.  Theoretically,  the  range  sidelobes  from  a  surface  scatterer  (and  its  first  periodic  image 
component)  can  also  be  reduced  by  decreasing  § f  beyond  the  required  value  of  5f  <c  /  2r  max ; 

however,  it  is  seen  that  this  approach  only  produces  results  with  limited  quality.  Better 
perfonnance  and  higher  algorithm  efficiency  are  realized  with  the  windowing  method. 

Also  in  order  are  a  few  notes  to  summarize  the  important  differences  between  the  single-  and 
multi-aperture  cases — in  terms  of  imaging  performance  as  conveyed  by  the  PSF  formulated.  As 
observed  with  the  set  of  parameters  applied  herein,  multi-aperture  sensing  leads  to  lower 
sidelobe  levels  in  the  cross-range  and  elevation  directions;  however,  the  range  sidelobe  levels 
and  range  resolution  are  comparable  to  those  of  the  single-aperture  case.  In  the  cross-range 
direction,  the  multi-aperture  configuration  gives  better  resolution  for  scatterers  situated  away 
from  the  broadside  direction  and  enables  attenuation  of  the  ambiguous  image  components.  (For 
the  multi-aperture  case,  the  ambiguous  component  for  each  target  now  manifests  as  a  series  of 
peaks — with  each  peak  corresponding  to  one  aperture.  Hence,  the  main-lobe-to-ambiguous- 
response  ratio  is  2 Nx .)  In  the  9  direction,  the  resolution  is  the  same  for  the  two  configurations. 

In  the  9'  direction,  better  resolution  is  noted  with  multi-aperture  imaging  for  scatterers 
positioned  at  azimuth  angles  away  from  (j)  =  ±n  /  2  ,  toward  broadside. 

2.2  Subsurface  Scatterer 

Equations  5  and  7  are  applicable  for  on-surface  targets.  Employing  the  appropriate  Green’s 
function  from  Reference  2  for  subsurface  observation  points,  the  PSF  for  a  shallow-buried  point 
target — for  the  multi-aperture  case — can  be  shown  to  be  of  the  same  form  as  Eq.  7,  except  the 
first  sin(-)  /  sin(-)  term  must  be  modified  as 
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and  ec  (fc)  is  the  complex  relative  dielectric  constant  of  the  ground  (at  the  center  frequency). 
Similarly,  the  first  sin(-)  /  sin(-)  term  of  Eq.  5  can  be  replaced  with  Eq.  12  sans  the 
Ax(sin6)cos0  -sin<9s  cos0s)  /  4  tenn  to  obtain  the  subsurface,  single-aperture  PSF.  Using  these 

subsurface  functionals,  for  sensing  and  imaging  at  the  grazing-angle  regime,  it  can  be  shown  that 
the  frequency  sampling  criterion  and  the  range  and  cross-range  resolutions  follow  the  same 
standard  forms  as  previously  established  for  the  on-surface  case,  while  the  resolution  in  6 
approximately  takes  the  fonn 
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where,  for  convenience,  A  (  /<  )  should  be  estimated  with  only  its  real  part — assuming  the 

imaginary  part  is  relatively  small.  A  closed-fonn  expression  for  the  resolution  in  6'  for  the 
subsurface  case,  however,  cannot  be  so  readily  obtained.  A  more  meaningful  numerical  analysis 
of  the  resolution  in  elevation  (or  in  depth)  can  also  be  reached  by  inspecting  the  behaviors  of  the 
exact  PSF  along  constant  phase  lines  below  the  ground. 


3.  Conclusions 


Analytical  closed-form  formulations  for  the  single-  and  multi-aperture  PSF  of  a  forward-looking 
radar  imaging  system  are  derived  in  this  work  for  both  on-surface  and  subsurface  point 
scatterers.  It  is  shown  that  the  multi-aperture  configuration  has  reduced  cross-range  sidelobe 
levels  and  better  cross-range  resolution  (away  from  the  broadside,  or  forward,  direction)  as 
compared  to  the  single-aperture  configuration,  and  the  ambiguous  image  component  in  the 


6 


azimuth  direction  can  be  suppressed  by  an  integration  of  the  multi-perspective  scattering  data. 
However,  the  sidelobe  levels  and  resolution  in  the  range  direction  are  the  same  for  both  the 
single-  and  multi-aperture  cases.  The  approximate  expressions  presented  herein  are  appropriate 
for  a  bistatic  sensing  configuration  in  which  the  excitation  is  supplied  by  the  two  transmitters 
located  at  the  ends  of  a  linear  array — a  sensing  geometry  that  is  consistent  with  the  one 
employed  by  an  existing  forward-looking  radar  setup;  nevertheless,  the  derivations  follow  a 
framework  that  can  be  applied  to  other  sensing  configurations  not  explicitly  treated  in  this  work. 
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